Effect of nuclear quadrupole interactions on the dynamics of two-level systems in 
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The standard tunneling model describes quite satisfactorily the thermal properties of amorphous 
solids at temperatures T < IK in terms of an ensemble of two-level systems possessing logarith- 
mically uniform distribution over their tunneling amplitudes and uniform distribution over their 
asymmetry energies. In particular, this distribution explains the observable logarithmic temperature 
dependence of the dielectric constant. Yet, experiments have shown that at ultralow temperatures 
T < 5mK such a temperature behavior breaks down and the dielectric constant becomes tempera- 
ture independent (plateau effect). In this letter we suggest an explanation of this behavior exploiting 
the effect of the nuclear quadrupole interaction on tunneling. We show that below a temperature 
corresponding to the characteristic energy of the nuclear quadrupole interaction the effective tunnel- 
ing amplitude is reduced by a small overlap factor of the nuclear quadrupole ground states in the left 
and right potential wells of the tunneling system. It is just this reduction that explains the plateau 
effect . We predict that the application of a sufficiently large magnetic field B > 10T should restore 
the logarithmic dependence because of the suppression of the nuclear quadrupole interaction. 

PACS numbers: 6143.Fs, 77.22.Ch,75.50Lk 
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Low temperature properties of amorphous solids have 
been associated with two-level systems (TLS) |l|,|2j. -Each 
TLS can be described by the tunneling amplitude Ao 
and the two-well asymmetry energy A. The tunneling 
amplitude Ao couples two energy minima and it is expo- 
nentially sensitive to external and internal parameters of 
the system. This results in the logarithmically uniform 
distribution!! 0,0,0 

P(A,A ) = -L^P = const. (1) 

Equation Q leads to a universal temperature and time 
dependence of various physical properties of amorphous 
solids, thus making "glassy" behavior easily recognizable. 
One of the signatures of the glassy state is a logarithmic 
temperature dependence of the dielectric constant and 
the sound velocity 0, . 

Recent experimental investigations on many amor- 
phous solids at very low temperature have revealed 
a number of qualitative deviations from the standard 
tunneling model. In particular, it was demonstrated 
H SH that for T < 5mK the dielectric constant 
is approximately temperature independent rather than 
logarithmically temperature dependent. These results 
conflict with the logarithmically uniform distribution Eq. 
JTJ of TLS over tunneling amplitudes. To resolve the 
problem one can assume that the distribution of TLS 
has a low energy cut-off at A . min ~ 1 — 3mK (See 0, 
Fig. P). Although this assumption fits the experimen- 
tal data quite good, it is in conflict with the observation 
of very long relaxation times in glasses 0]. These times 
(a week or longer) require much smaller tunneling am- 
plitudes than lmK (remember that the TLS relaxation 



time is inversely proportional to its squared tunneling 
amplitude) . 

We suggest an explanation for this controversial behav- 
ior using recent work 0, 0, 0| , where it was proposed 
that the nuclear quadrupole interaction affects the prop- 
erties of TLS's at low temperatures due to a mismatch of 
the nuclear quadrupole states in the two potential wells 
. The case of higher temperatures T > 5mK has been 
considered in [ill fl3 | . The significance of the nuclear 
quadrupole interaction has been proven experimentally 
in glycerol glass 0] . This interaction also helps to under- 
stand the unusual magnetic field dependence of dielectric 
properties in non-magnetic dielectric glasses 0, 0] • 

The nuclear quadrupole interaction splits the levels of 
a well and the effective tunneling between pairs of them 
is Ao* < Ao, i.e., gets smaller than the initial value Ao. 
We demonstrate that this reduction is of decisive impor- 
tance for TLS with tunneling amplitude Ao smaller than 
the characteristic nuclear quadrupole interaction energy 
A*. The contribution of these TLS's to the dielectric 
constant is substantially suppressed. This can explain 
the experimental observations. 

First, we review the key data concerning the contri- 
bution of tunneling systems to the resonant part of the 
dielectric susceptibility for T < 50mK . This con- 
tribution can be described as the adiabatic response of 
the quantum system in the ground state to a slow ex- 
ternal field. The application of an external electric field 
to the tunneling system in the ground state produces a 
polarization along the field. For instance, in the two- 
level model the external field E affects the TLS Hamil- 
tonian through the change in the asymmetry energy, i.e., 
-A s x - (A - /xE)s z [! 0. The ground state has the 
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energy E g = — (1/2)^ A 2 + A§, which results in the sus- 
ceptibility Xg = -^^Eg/dA 2 = i r 2(A .+|2 ) 3/2 - 

A contribution of excited states to the susceptibility 
can be negative. In a two-level system the contribution of 
the excited state with the energy E e = (l/2)yj A 2 + Aq 
is — Xg- When the temperature is less than the minimum 
energy splitting between the ground and excited state 
E m i n (Ao) — Ao, the excited state contribution can be ne- 
glected. Then, the contribution of the TLS's with given 
Aq > T and different A to the susceptibility is 
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where W 3> Ao, T is the maximum TLS asymmetry. For 
large enough A one has E g (A) ~ — | A | /2 and the inte- 
gral 10) equals unity. This result is valid for the ground 
state of the tunneling system irrespective of the presence 
of the nuclear quadrupole interaction. Then, one should 
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make averaging of -P%- over Ao using the distribution Eq. 
P|l with the constraint that the minimum energy split- 
ting E m in between the ground- and first excited state 
exceeds T. This is realized when Ao > T. Otherwise, 
the positive dielectric susceptibility of the ground state 
is compensated by the negative contribution of the ex- 
cited state and the total contribution can be neglected. 
Then one may express the TLS resonant susceptibility to 
logarithmic accuracy in the form 
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e(£ mi „(A )-T), (3) 



where Q(x) is the Heaviside step function. In the stan- 
dard TLS model E min (A ) = A 0, Q and we obtain 
the well-known result xtls ~ (P/j? /3)ln(W/T). Equa- 
tion is still applicable in the low temperature limit, 
where the nuclear quadrupole interaction is significant. 

To describe the effect of the nuclear quadrupole inter- 
action on the dielectric constant, we consider the change 
in TLS ground state properties induced by this interac- 
tion. The adiabatic (resonant) dielectric susceptibility 
X oc A~ 3 (see Refs. H 0) is due to the TLS's with 
small asymmetry energies A < Ao. Moreover, as dis- 
cussed previously, the minimum energy splitting between 
the ground and first excited state takes place at A = 0. 
Therefore, in what follows we study this particular case. 

Consider a two-level system formed by n tunneling 
atoms possessing a nuclear spin I > 1 and consequently 
a nuclear quadrupole moment. The tunneling system 
Hamiltonian can be expressed through the standard TLS 
pseudospin Hamiltonian — Aos x and the quadrupole in- 
teractions, which are H r and Hi in the left well (s z — 
— 1/2) and the right well (s z = 1/2), respectively. The 
mismatch of the nuclear quadrupole interactions between 
the two wells makes the TLS properties sensitive to an 
external magnetic field 10]. We also assume that the 



choice of A = implies identical ground state ener- 
gies of nuclear spins in the right and left wells, i.e., 
< gl | H t | gl >=< gr | H r | gr >= E l ° u c adr , where 
\gl > (\gr >) is the ground nuclear spin state in the left 
(right) well. The modified TLS Hamiltonian takes the 
form 



h TLS = -A s x + {H r + H t )/2 + (H r - H t )s* 



(4) 



Two different realizations of the ground state are pos- 
sible for the Hamiltonian Eq. depending on the rel- 
ative size of the TLS tunneling amplitude Ao and the 
nuclear quadrupole interaction energy A*. If Ao > A*, 
then the hybridization is strong and the tunneling par- 
ticle is equally shared between both wells and one can 
set s x = 1/2, s z pa in Eq. Q. Then, the ground 
state energy is estimated as Eg yb = — A /2 + E ^ adr , 
where E q ^ adr is the ground state energy of the nuclear 

spin in the hybridized Hamiltonian (H r + iJ;)/2. If 
Ao < A*, the ground state of the system corresponds 
to either the localized left well ground state or the right 
one E l ™ adr = <gl\Hi\ gl >=< gr | H r | gr >. The 
tunneling term in Eq. Q leads to a small correction to 
this energy given by 



SE = -t^Ao/2, ry* =< gl\gr > 



(5) 



The parameter 77* < 1 characterizes the average overlap 
integral of the nuclear spin ground states in the left and 
right wells. We assume that ?y„ -C 1. For this reason, cor- 
rection (jSJ) is ignored below when estimating the energy. 

The energy splitting between the ground and the first 
excited state, formed by the symmetric and antisymmet- 
ric superpositions of the localized ground states in the 
left and right wells is given by the modified effective tun- 
neling amplitude 



A * = r;*A = E min (A ), 



(6) 



which also defines the minimum energy splitting of the 
ground and first excited states i? m j n (Ao), because other 
excited states can be neglected (see below Eq. (JSJ and 
the related discussion). 

We have found the two possible realizations for the 
ground state of TLS's affected by the nuclear quadrupole 
interaction, a hybridized and a localized one. For the 
parameter A* introduced above one has 
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In the case Ao > A*, i.e., the TLS is hybridized between 
the left and right wells and the nuclear quadrupole in- 
teraction can be approximately ignored. In the opposite 
case Ao < A*, the nuclear quadrupole interaction local- 
izes the tunneling particle in one of the wells leading to 
a strong reduction of the effective tunneling amplitude 
and replacing Aq by Aq* = 77* Aq Aq. The nuclear 
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spin ground state in each well can be treated as non- 
degenerate ^f|. Then, one can approximately define the 
ground state and the lowest excited state as superpo- 
sitions of these two lowest energy states. The ground 
state of the nuclear spin in each of the two wells can 
be considered as non-degenerate and one may construct 
the ground state and the lowest excited state as a super- 
position of unperturbed ground states in the two wells. 
Higher excited state contributions are neglected. This is 
justified by calculating their contribution to lowest order 
perturbation theory. We find a correction factor c to the 
ground state amplitude as 
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We determine the parameter regime for which the second 
term can be neglected as follows. Assume that each TLS 
contains n simultaneously tunneling atoms. The lowest 
excited states in each well are separated from the ground 
state by some characteristic energy b ~ huj$ where loq is a 
frequency of the nuclear quadrupole resonance. The next 
group of states is separated by the energy interval ~ 2b. 
The vast majority of the states have energies exceeding 
that for the gro und state by the energy nb ~ A*/2 (see 
Eq. J7J and [l6j). Because of the large statistical weight, 
these states provide the main contribution to Eq. JHJ. 
So, we may replace the denominator in Eq. JHJ by A*. 
Then, the sum of the |< gl\i >| 2 in the numerator can 
be approximated by 1. The requirement that the second 
term in Eq. (JSJ is small results in the condition 



A < A, 



(9) 



Thus c ss 1 when A < A* so in that case our considera- 
tion can be restricted to the two lowest energy levels. 

Consider the adiabatic dielectric susceptibility at dif- 
ferent temperatures. For high temperature T > A* the 
quadrupole interaction is not significant and we may use 
the result of the standard tunneling model 

Xtls(T) « (P» 2 /3)ln(W/T), K<T<W. (10) 

Now consider the limit of low temperatures, where the 
main contribution to the temperature dependence comes 
from the TLS's with modified tunneling splitting Eq. 
According to Eq. J3J, the TLS contribution to the dielec- 
tric constant can be estimated by the logarithmic inte- 
gral with the lower limit Ao, mim defined by the condition 
E m i n (Ao, m in) == T. At low temperature this minimum 
energy splitting is given by the renormalized tunneling 
amplitude Ao* = 77* Ao (see Eq. {JJ) and the integral in 
Eq. (O should be cut-off at Ao jm i„ = T/77* rather than 
at Ao jm m = T as in the standard tunneling model. This 
leads to a negative correction to Eq. (|10|l of the form 



The condition for the applicability of Eq. is given by 
the criterion Eq. ((3) which has to be valid for tunneling 
system with the smallest involved tunneling amplitudes, 
in particular Ao )ITO *„ < A*. Using Ao )Tn j„ < T/77* we 
obtain the criterion for the validity of Eq. Ijllf) 



T < 77* A* 



(12) 



The value of Eq. Qlljl at the highest admissible temper- 
ature 77* A* coincides with the value of Eq. l|l(Jf) a t the 
lowest admissible temperature A*. Since the adiabatic 
dielectric constant cannot decrease with decreasing tem- 
perature, because excited states have a negative effect 
on it, we expects that it does not change at intermediate 
temperatures (plateau regime) 

XTLs{T)~{P^/3)ln(W/K), A*?7*<T<A*. (13) 

The plateau can also be interpreted as the consequence 
of the pseudogap in the distribution of effective tunneling 
amplitudes Ao* within the range A* 77* < Ao* < A* |16| . 

To compare the theory with the experiment let us es- 
timate the parameter A*. We can use the approximate 
expression obtained above, i.e., A* ~ nb where n is the 
number of atoms per two level system and the energy b 
has been defined before. The majority of glasses show- 
ing the plateau in the temperature dependent dielectric 
constant 0, 0, contains K, Na, Al, Ba atoms which 
possess a relatively large nuclear quadrupole moment (see 
13, flfij]). The nuclear quadrupole splitting energy b mea- 
sured for these nuclei in different environments varies be- 
tween QAmK and lmK 16]. To obtain A* ~ 2nb ~ 5mK 
(see Eq. (00 one can suppose that n « 4 atoms are in- 
volved in the tunneling motion. This estimate looks quite 
reasonable taking into account the strong effect of each 
tunneling atom on its neighboring atoms. Tunneling of 
the single atom shifts equilibrium positions of its neigh- 
bors. During the tunneling event the neighbors move to 
new equilibrium positions simultaneously with the tun- 
neling atom. Therefore they participate in tunneling and 
should be included into the definition of the parameter 
n, which can thus become larger than 1. The overlap 
parameter 77* =< lg\rg > is taken between the ground 
states of n non-interacting nuclear spins and estimated 
as 



77* = 7/ 



(14) 



Xtls(T) « (Pn 2 /3)(ln(W/T) - ln{l/n*). 



(11) 



where 77 is the ground state overlap integral for one spin. 

For a crude estimate of 77 w e consider the orientational 
glass (KBr)i- x (KCN) x [17| where the CN group ro- 
tates between different equilibrium positions by an an- 
gle (j> = cos _1 (l/3) w 0.123 3]. For the sake of sim- 
plicity consider the simplified model of a nuclear spin 
1 = 1 with a quadrupole interaction described by the ax- 
ially symmetric Hamiltonian [13 6(J2-2/3) andfe> 0. 
The index a denotes the nuclear quadrupole axis. Then 
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-nB=0mK, B=0 

- modeling of experiment 
' B=0, A„ . =1.8mK[5] 

Omin L J 

-nB=2mK, B-10T 
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FIG. 1: Temperature and field dependence of TLS dielec- 
tric constant. The solid line shows the temperature depen- 
dence of the dielectric constant, calculated using parameters 
A* = 4mK, r\ — 0.33, n = 4 and the tunneling amplitude Ao 



replaced with Ao* = An?7 
rived in Ref . HF 



n(l-(A /A„) 2 ) 3 / 2 ) 



at Ao < A* as de- 



The artificially generated solid line with cir- 
cles represents experimental data [3, which can be described 
within the standard tunneling model with the tunneling am- 
plitude cutoff Aomin = 1.8mK (representative value, see ||). 
It is clear that our approach does not deviate from the ex- 
periment in the experimental temperature range T > lmK. 
Dielectric constant in the external magnetic field (dot-dashed 
line) is found using the numerical modeling for axially sym- 
metric nuclear quadrupole interaction Ht^ . 



one may estimate a single nuclear overlap integral as 
r\ w cos((f>) « 0.33. If a TLS consists of n = 4 tunneling 
atoms, the effective overlap integral for two ground states 
can be estimated as tj* ss 10~~ 2 . In this case 77* A* « 50[iK 
and the plateau should show up in the temperature range 
from T down ~ 50fiK to T up - AmK (see Eq. (|T3J)). For 
T < 50fiK the logarithmic temperature dependence of 
the dielectric constant should be restored and resemble 
the one at high temperatures (see Fig. Note that 
the measurements of Tdown and T up ) will permit us to 
estimate the number of atoms participating in TLS as 
n ~ ln(T up /Tdown)- This is important for understanding 
the microscopic nature of tunneling systems. 

The theory can be verified by the measurements of the 
dielectric constant at low temperatures T < A* under the 
strong magnetic field, affecting the nuclear magnetic mo- 
ments. When the field strength becomes comparable to 
the nuclear quadrupole interaction, it changes the spin 
quantization axes in both left and right wells aligning 
the nuclear spins along the field. Accordingly the mis- 
match of spin quantization axes is reduced by the mag- 
netic field and the overlap integral ry* increases and ap- 
proaches unity in the large field limit. Thus we predict 
that a very strong magnetic field eliminates the influ- 
ence of the nuclear quadrupole interaction on the dielec- 
tric constant and restore the standard tunneling model 
behavior. The critical magnetic field strength at which 



the nuclear quadrupole effects are eliminated can be esti- 
mated by requiring that the Zeeman splitting in the field 
is comparable with the energy A* ~ AmK . This expecta- 
tion agrees with our numerical estimate |l6| . Assuming 
that for the typical nuclear magnetic moment /1 one has 
/j,B ~ lmK at B ~ 5T, one can estimate the magnetic 
field necessary to eliminate the plateau in the dielectric 
constant to be B ~ 10T. The calculations of the temper- 
ature dependence of the dielectric constant at zero and 
large B [16( are shown in Fig. Q] together with the plot, 
representing the experimental data |5(. 

In summary, we explained the plateau in the temper- 
ature dependence of the dielectric constant by the local- 
ization of the tunnelling particle by nuclear quadrupole 
interactions. We predict that an external magnetic field 
B > 10 T will restore the logarithmic temperature depen- 
dence of the dielectric constant in the whole temperature 
range. Such an experiment will verify the theory. 
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supported by the program of Russian Scientific school. 



[1] 

[2] 
[3] 

[4] 
[5] 

[6] 



[7] 



[8] 
[9] 

[10] 

[11] 
[12] 

[13] 



[14] 
[15] 



[16] 



P.W. Anderson, B.I. Halperin, CM. Varma, Philos. Mag. 
25, 1 (1972) 

W.A.Phillips, J.Low Temp.Phys. 7, 351 (1972). 

S. Hunklinger, A. K. Raychaudchary, Progr. Low Temp. 

Phys. 9, 267 (1986). 

W. A. Phillips, Rep. Prog. Phys. 50, 1567 (1987). 

S. Rogge, D. Natelson, B. Tigner, D. D. Osheroff, Phys. 

Rev. B 55, 11256 (1997). 

A. L. Burin, D. Natelson, D. D. Osheroff, Yu. Kagan, 
in "Tunneling Systems in Amorphous and Crystalline 
Solids", ed. P. Esquinazi, p. 243 ( Springer, 1998) 
P. Strehlow, W. Wolhardt, A. G. M. Jansen, R. Haueisen, 
G. Weiss, C. Enss, S. Hunklinger, Phys. Rev. Lett. 84, 
1938 (2000). 

C. Enss, Physica B 316, p. 12 (2002). 

D. J. Salvino, S. Rogge, B. Tigner, D. D. Osheroff, Phys. 
Rev. Lett. 73, 268 (1994). 

A. Wiirger, A. Fleischmann, C. Enss, Phys. Rev. lett. 89, 
237601 (2002). 

D. Bodea, A. Wiirger, J. Low Temp. Phys. 136, 59 (2004) 
I. Ya. Polishchuk, P. Fulde, A.L. Burin, Yu. Sereda, D. 
Balamurugan, J. Low Temp. Phys. 140, 355 (2005) 
To our knowledge the only one glass shows no plateau 
in the dielectric constant, which is mylar, containing no 
nuclei with spin I > 1/2 Q. This experimental result and 
similar results in other systems support our theory |lq|. 
P. Nagel, A. Fleischmann, S. Hunklinger, C. Enss, Phys. 
Rev. Lett. 92, 245511 (2004). 

Nuclear spin Hamiltonians for n-atomic tunneling sys- 
tem can be represented as the sum of n single atom 
quadrupole interactions [iHf . For integer nuclear spin the 
ground state is non-degenerate. The same is true for the 
whole Hamiltonian. If the nuclear spin is half-integer, it 
has a double-degenerate eigenstate and the ground state 
of n spins is 2™ times degenerate. One can show [l6j that 
these ground states can be separated into independent 
2 n pairs, which can be treated as a non-degenerate ones. 
A. L. Burin, I. Ya Polishchuk, P. Fulde, Yu. Sereda, sub- 
mitted to Phys. Rev. B; cond-mat/0509810 



■5 



[17] J. N. Dobbs, M. C. Foote, A. C. Anderson, Phys. Rev. B 
33, 4178 (1986). 



